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Abstract 

Spin transport properties of a spin-polarized two-dimensional electron gas 
are studied in the presence of electron-electron interactions. Longitudinal 
and transverse spin diffusion coefficients are calculated with the quantum 
transport equation which includes many-particle effects in the random phase 
approximation. We find that the e-e scattering, which does not contribute to 
charge drift mobility, has a significant contribution to the spin diffusion. Thus, 
e-e interaction causes non-negligible effects on the operations of proposed 
devices dependent on spin transport in semiconductor heterostructures. 
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I. INTRODUCTION 



Spin dependent electronic transport phenomena in semiconductor nanostructures are ex- 
tensively studied in recent years, with motivations to develop novel electronic and optical 
devices based on carrier spin dynamics such as spin transistorsB, polarization dependent 
optical modulators, and even quantum computers. Recent experimental progress in trans- 
porting electronic spins! over the distance of 100 /zm and injecting spinsll from ferromagnetic 
materials to compound semiconductors has been accelerating the prospect of these devices. 

Charge transport is extensively studied to understand the operations of devices. In 
the conventional devices where the electron distribution is spatially homogeneous and the 
electrons move in the same direction irrespective of spin directions, the charge transport is 
expressed in terms of drift mobilities. When an electronic band is parabolic and umklapp 
processes are negligible, electron-electron (e-e) scatterings do not directly affect mobilities 
since the e-e scattering conserves the total momentum of the system.! It influences the charge 
transport only indirectly by modifying momentum distribution functions and screening other 
scattering processes of electrons. Mobility is determined by e-phonon and e-ionized impurity 
scatterings. In high-quality, modulation-doped heterostructures of III-V semiconductors, 
electron mobility increases as temperature decreases and reaches the maximum value limited 
by the scattering with remote dopantsJl 

In contrast to the charge transport, spin-up electrons move against spin-down electrons 
or vice versa in the spin transport. Thus, in addition to e-ionized impurity scattering, we 
expect that the Coulomb interaction between spin-up and -down electrons also contributes 
to the spin transport.! (The Coulomb interaction between two electrons with the same spin 
direction does not directly influence the spin transport as it conserves the spin current.) As 
we show later, the e-e interaction is as important as the e-ionized impurity scattering to the 
spin transport. 

The spin transport is also influenced by the many-particle effects of interacting electrons. 
An electronic band dispersion is renormalized, giving the effective mass of an electron some- 
what different from the band effective mass, and the bare Coulomb interaction is screened, 
leading to the e-e interaction with a finite range. These many-particle effects affect collision 
dynamics and modify transport coefficients. 

The many-particle effects also have consequences on spin dynamics itself. The spin axis 
of a noninteracting particle precesses around external magnetic field B e:r . In the interacting 
system, each electron feels effective magnetic field B mo /(r) due to the molecular field of 
surrounding electrons. The effective field is related (almost proportional) to the local spin 
density and parallel to its direction as explicitly shown in our analysis below. Thus in a 
spin polarized electron gas, a spin axis precesses around B ex + B mo ;. This induces a "spin- 
rotation term" in the equation of motion for the spin current, leading to a anomalous spin 
diffusion coefficients and spin wavesll'l which were experimentally observed in 3 He, 3 He- 4 He 
mixture, and in atomic HydrogenBEII 

In this paper we present a theory and numerical calculations of spin diffusion coefficients 
in a two-dimensional electron gas (2DEG) interacting through the Coulomb force. In our 
previous publication^ we adopted the Hartree-Fock approximation (HFA) to incorporate the 
many-particle effects in the lowest order. Here we adopt the random phase approximation 
(RPA) to improve the treatment on many-particle aspects of electrons and obtain more 
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reliable numerical results. We compare the contributions from e-e scattering and e-ionized 
impurity scattering to the spin diffusion in semiconductor heterostructures and discuss their 
relative importance. 

This paper is organized as follows: In Sec. II, our calculation scheme is briefly described. 
In Sec. Ill the numerical results for the spin transport coefficients are shown with detailed 
discussions. Conclusions are given in Sec. IV. The differences in mathematical expressions 
from our previous calculation are summarized in Appendix. 



II. THEORY 

Our theoretical model is already presented in our previous publication^ and we describe 
it only briefly here. Several relevant expressions arising from the RPA will be given in 
the appendix. We consider two-dimensional electrons interacting through the Coulomb 
interactions. Other interactions such as electron-phonon, electron-ionized impurity or the 
exchange interaction of electrons with magnetic impurities are not included in the present 
study. We assume that electronic spins are polarized in the growth direction, z-axis (normal 
to the 2D plane), in the quasi-thermal equilibrium, and a direction of local spin polarization 
can be slightly tipped away from z-axis by applying r.f. field or by external spin injection. 
The total sheet density of spin-up and spin-down electrons, n(r,t) = n+(r, t) +n_(r, t), is 
assumed to be spatially and temporally constant,and the spin density, M(r, t) = n + (r,t) — 
n_(r,t), varies slowly on a hydrodynamic scale.Ea The spin polarization is given by P = 
(n + — n_)/ (n + + nJ). Spin relaxation is not included in the present analysis. 

A moderate external magnetic field, B ex , can be applied perpendicular to the 2D plane. 
We assume that the field is so weak that it only affects spinor space, leaving electron orbital 
motion intact. Otherwise the orbital motion is quantized into the Landau levels, and the 
system turns into the quantized Hall state, to which our present theory does not apply. The 
estimation in our earlier work0 shows that the weak field condition is satisfied when B ea , < 1 
T. When this condition is satisfied, the external field can trivially be erased from the theory 
by using a rotating frame. 

The quantum transport equation for the electron distribution functions n k is given byEl'i 



dn p (r,t) 1 j de p (r, t) dn p (r,t) \ _ 1 ( de p (r, t) dn p {r,t) \ %_ 
dt 2 1 dp ' dr J 2 1 dr ' dp J % 



£p(r,t),n p (r,t)] 
/ dn p (r,t)\ 



Of 



(2.1) 



col. 



where the electron energy e k and the distribution function n k are 2x2 matrices in spinor 
space, denoted by underlining. The diagonal components of n k represent the population 
of spin-up and spin-down electrons and its off-diagonal elements represent the couplings 
between spin-up and spin-down states. The transport equation is similar to the Boltzmann 
equation except that there are anticommutators in the drift terms (the second and the third 
term on the left) and a commutator in the last term on the left. This term in the commutator 
is called the "spin-rotation term" , representing the effect of the external and molecular fields. 

The electron (quasiparticle) energy e k is composed of a kinetic energy term, a coupling 
to the external field and a self-energy caused by the presence of other electrons. In our 



3 



earlier work,0 we adopted the HFA to calculate an electron self-energy since the simple 
analytic expressions of the HFA considerably reduces the load of numerical calculations. 
But the HFA does not describe the many-particle properties of 2DEG appropriately: The 
energy dispersion of a quasiparticle is not given properly with a vanishing quasiparticle 
effective mass at Fermi surface due to the infinite range of the bare Coulomb interaction. 
Since the effective mass is one of the important parameters in our theory, a certain degree 
of arbitrariness exists in our previous calculations. And since the HFA does not include 
correlation terms, the self-energy for a spin-up electron does not contain the contribution 
from spin-down electrons. This leads us to an unrealistic situation when the spin polarization 
is high, i.e., provided that we prepared 2DEG with all spins pointing downward, and we put a 
spin-up electron into this 2DEG, then the energy dispersion of this spin-up electron remains 
that of a noninteracting electron in spite of the presence of many spin-down electrons. This 
prevented us from calculating the diffusion coefficients of highly polarized 2DEG in our 
earlier work. 

In the present study we improve our calculation by adopting the RPA. It takes into 
account exchange and correlation (though not completely), and we can correct the problems 
of the HFA mentioned above. The RPA self-energy for an electron gas is discussed in many 
text books (e.g., see Ref. [TJ] and |T5p, and its detailed expressions for a two-dimensional, 



multi-band electron system are given by VinteiE§ and DasSarma et. We use the stan- 
dard RPA self-energy without vertex corrections. The Matsubara function (imaginary time 
Green's function) for the self-energy, which is a 2x2 matrix in spinor space, is written as 

h^ A (k,ik m ) = -^^^^(- q ,-^ n )G Qi/3 (k + q,^ m + ^) ; (2-2) 

where a and (3 are indices for spin, G° a f3 is the noninteracting electron propagator. (G° a * is 
finite only when a = (3.) The RPA screened Coulomb interaction V^ A satisfies the Dyson 
equation, 

V3f A (q,«&0 = Va,p(<l) +Y.V<*MKM,nn)V s J A {cL,iqn), (2.3) 

8,5' 

where V a ^(q) is the bare Coulomb interaction. The lowest order polarization (bubble dia- 
gram) is finite only when 5 = 5', given by 

Ps,s(^ i( ln) = E E G h(P, iPm)G°sAP + q, iPm + iqn). (2.4) 

P P iPm 

The real-time functions are obtained by the analytic continuation after the momentum 
integration and the summation over Matsubara frequencies. The quasiparticle energy dis- 
persion is obtained from the real part of the self-energy following the prescription suggested 
by Rice.0 

In Fig. 1 we show energy-momentum dispersions of spin-up and -down electrons in 
2DEG with spin-polarizations P = 0.1 - 0.9 and with total electron sheet density at 2 
x 10 11 cm ~ 2 . The electron energy is measured relative to the zero of the noninteracting 
electron dispersion. One of the important many-particle effects in the interacting electrons 



4 



is the reduction of one-electron energy (Bandgap renormalization) . Since n + > n_, the 
renormalization is larger in the spin-up electron. The energy differences between the spin- 
up and the spin-down electron increase with spin polarizations. When P < 0.5, the energy 
dispersions of both spin-up and -down electrons are well approximated by parabolic band 
dispersions. We obtained the quasiparticle effective mass m* between 1.11m and 1.19m 
depending on spin polarizations, (m is a band effective mass. We used m = 0.067mo for 
the conduction band of GaAs.) In numerical computations of transport coefficients, we have 
used the parabolic band approximation for the energy dispersions. 

The quantum transport equation eq.( |2.1p is solved by the Chapman- Enskog expansion 
to the lowest order using variational solutions.0 The final expression of the equation for the 
spin current is obtained by multiplying the transport equation by the particle velocity and 
integrating over the momentum as,EJ 

d J (r t) dM(r t) de 

Q t - ! + e a 1 A (k F ±) + M(r, t)-^-B (k F± ) + HsrM x J CTji + lgy B ex x J v%i 

= "tr\\,i "tr±,ii (2-5) 

where .j is a spin current representing the flow of the spin component a in the spatial 
direction j = x, y, z, e a is a unit vector in the direction of spin polarization, and M = Me a . 
The factor fisR, which is a function of k F ±, represents the contribution of the molecular 
field, and 7 TO = —gfis/h is a gyromagnetic ratio. The coefficients A (k F ±), B (k F ±) and fisR 
will be given in the appendix. 



III. RESULTS AND DISCUSSION 
A. Spin diffusion coefficients 

We show the results of numerical calculations for spin diffusion coefficients. Since a spin 
current is a tensor defined by the two vectors; the direction of spatial flow and the direction of 
spin, we should distinguish longitudinal (Dm) and transverse (D±) spin diffusion coefficients. 
(See Fig. 2) In the longitudinal spin diffusion, consider the system where the direction of 
spin polarization is aligned in z everywhere but its magnitude has gradient in the direction 
x. Then the longitudinal spin current flows in the spatial direction x whose spin direction 
is in z, parallel to the spin polarization. While in the transverse spin diffusion, consider the 
system where the magnitude of spin polarization is the same everywhere but its direction is 
tilted gradually from z to x by a small angle when we go from x = Xq to Xq + dx. In this 
case the transverse spin current flows in the spatial direction x with the spin direction in x, 
perpendicular to the spinoolarization. These two diffusions show different magnitudes and 
temperature dependence.!^ 

The temperature dependence of D\\ and is plotted in Fig. 3 and 4, respectively, with 
0.01 < P < 0.5. We calculate in the degenerate region, T < 20 K. The Fermi temperature is 
83 K in the unpolarized 2DEG with the sheet density 2 x 10 11 cm~ 2 (Fermi wave number k F 
= 0.112 nm" 1 ). The magnitudes of the diffusion coefficients are somewhat larger than those 



calculated previously in Ref. O but their temperature dependencies are identical. We believe 
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that the present results calculated with the RPA should be more accurate in predicting 
the transport coefficients. As shown in Fig. 3, Dn increases at low temperatures and 
diverges as T~ 2 when T approaches zero. A close inspection shows it deviates slightly from 
T~ 2 dependence (bulging downward). Actually we find D\\ oc {E F / ksT) 2 /\a.{E F / k B T) as 
expected of the "normal" transport coefficient in 2D degenerate Fermions. (The logarithmic 
correction factor is characteristic of a 2D system, i) While in Fig. 4, D± increases in the 
higher temperatures in a similar manner as D\\, but it departs from T~ 2 dependence and 
becomes constant as T — > at the larger P. 

The different temperature dependence of D» and D± is attributed to the phase space 
available in the collisions of two diffusion processes. In the longitudinal spin diffusion, a 
close inspection of the associated collision term shows that electrons in the vicinity of the 
Fermi surface within the width k B T only are allowed to participate in the collision due to 
the kinematics and the Pauli exclusion. Thus as the temperature is lowered, the number of 
electrons participate in the collision decreases, leading to the reduction of a collision rate. 
While in the transverse case, all electrons between the two Fermi surfaces of the spin-up 
and -down electrons participate in the collision. Thus the larger the spin polarization is, 
the more electrons are involved in the collision, leading to the finite and constant diffusion 
coefficients at T — * 0. This anomalous temperature dependence is observed experimentally 
in 3 He- 4 He system.0 



B. Electron-impurity scattering 

In the previous subsection we have discussed the spin diffusion of 2DEG in the presence 
of e-e scattering alone. The spin diffusion in actual samples of heterostructures and quantum 
wells (QWs) should also be affected by other scattering processes. In heterostructures of III- 
V semiconductors, a low temperature electron mobility, /i e , is dominated by the scattering 
with remote ionized impurities (dopants). We expect that the e-impurity scattering will also 
strongly influence the spin transport. We compare the spin diffusion coefficients limited by 
e-e scattering which we have calculated above with those limited by e-impurity scattering, 
and show that the e-e scattering bears substantial contribution to the spin transport in 
actual sample structures, although the e-e scattering has only an indirect consequence to 
the charge transport. We do not give the explicit calculations of the spin diffusion including 
the e-impurity scattering but we estimate its magnitude from experimental carrier mobilities. 

The highest mobility reported so far in GaAs/AlGaAs heterostructures exceeds 1 x 
10 7 cm 2 V -1 s _1 below 1 K.@ The charge diffusion coefficient D e arising from e-impurity 
scattering is related to the mobility /i e through the generalized Einstein's relation.^ For a 
degenerate system, we have 

pD 

= (1 + g-zVfcBT) log(1 + e „/k B T, 

fx e k B T 

where /i is a chemical potential. In the high quality samples cited above, the charge diffusion 
coefficient D e is estimated to be 9 x 10 4 cm 2 s _1 below 10 K. Assuming that the magnitude 
of the spin diffusion limited by e-impurity scattering is similar to the charge diffusion, we 
obtain the estimated value of 9 x 10 4 cm 2 s _1 for the spin diffusion, which is, compared 
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with Fig. 3, close to the spin diffusion limited by e-e scattering at the lowest temperatures. 
Thus the e-e scattering has a significant contribution to the spin diffusion, comparable to 
the e-impurity scattering. In moderate-quality samples with \i < 1 x 10 6 cm 2 V -1 s _1 , 
the relative weight of e-ionized impurity scattering should become larger than that of e-e 
scattering. 

C. Effective magnetic field 

In a many-electron system an individual electron interacts with surrounding electrons 
through the Coulomb interaction. This effect is represented in the Fermi liquid theory that 
an electron is subject to a molecular (mean) field produced by surrounding electrons. Due 
to the exchange part of the interaction, the energy of an individual electron depends on 
the spin directions of itself and neighboring electrons. Thus the molecular field acts as an 
effective magnetic field B mo 2 to the individual electron, and its spin precesses around B moZ . 
The consequence of B mo ; was first pointed out by Leggett and RicJS that the effective 
transverse spin diffusion coefficients measured by spin echo experiments should depend on 
the spin-tipping angle by the initial r.f. pulse (Leggett-Rice effect). The anomalous behavior 
of the transverse spin diffusion and the spin waves discussed in our previous publication!^ 
have its physical origin in the precession of flowing spins around the effective magnetic field. 
Consequences of the effective magnetic field can be seen explicitly in the last two terms on 
the left-hand side of eq. ( j2.5|) , /is^M x J CT j + 7 TO B ea; x J^ j. The second term expresses 
the spin precession around the external magnetic field, while the first term represents the 
spin precession around the effective magnetic field. The effective field is parallel to spin 
density M, and its magnitude depends on fi SR and M. We should notice that the precession 
occurs not on the spin density but on the spin current J CT , i.e., the precession occurs on an 
electron flowing into the region with its spin direction slightly tipped away from the local 
spin density. (If the direction of M is constant throughout the system, the effective field 
B mo ; has no net effect.) In Fig. 5 we plot the magnitude of the effective field in terms of spin 
splitting energy, AE mo i = Ti^srM [eV], (or the Larmor precession frequency) as a function 
of spin polarization P = M/n. The spin splitting shows a sublinear dependence on the spin 
polarization since [asr decreases with P. 

IV. CONCLUSIONS 

We have calculated spin diffusion coefficients of 2DES interacting through the Coulomb 
force. The improved approximation in the RPA leads to a more precise treatment of many- 
particle aspects of the electron properties compared with our earlier calculations using the 
HFA. The magnitudes of both the longitudinal and the transverse spin diffusion coefficients 
are larger than our previous calculations but their temperature dependencies are qualita- 
tively similar to the previous ones. 

Charge transport in heterostructures, at low temperatures, is dominated by the electron- 
ionized impurity scattering, while the e-e scattering has little consequences. In contrast, as 
we have shown, the spin transport is largely affected by the e-e scattering in high-quality 
samples. Thus we expect that the spin transport properties in semiconductor devices should 
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depend on spin-up and -down carrier distributions both in real and momentum space, in 
addition to the distributions of remote dopants. This fact should be taken into account in 
designing recently proposed spintronics devices. 
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APPENDIX A: EXPRESSIONS IN RPA 

The equations are shown in detail in our previous publicationO calculated in the HFA. 
In the present study, the mathematical expressions can be obtained by substituting the self- 
energy terms calculated in the RPA for those in the HFA. In contrast to the HFA self-energy, 
which gives an analytic expression, the RPA expressions require numerical integrations. We 
give the relevant expressions used in the present study below. 

The quasiparticle energy £°(r, t) in the local equilibrium (Eq. (2.4) in Ref. [11]) is obtained 
from the real part of the RPA self-energy following the prescription suggested by Rice,@ 

e° p± (r, t) = e Kin + Refr£^ PA (p, Two = e Kin ; r, t), (Al) 

where + (— ) for (1,1) ((2,2)) component, respectively, and the off-diagonal elements are 
zero. Sxin — p 2 /2m is the noninteracting energy of an electron, and h'E± PA (p, hu; r, t) is 
obtained from the imaginary-time function in Eq. (|2.2|) . In the equation of motion for spin 
current Ja-j, the coefficient A is a function of spin-up and spin-down electron densities (or 
the Fermi wave numbers) and is given by 

A(k rt)- 1 / VM) , n-(r,t) \ 

A(k F± ,r,t) - _ I + l (A2) 



where G± = G± — a± with 



m ± UE P =hk F± ,± OE P =hk F ±,± 



The coefficient B is 



27ih 2 ' a± dn + dn 



7ih 2 n . 1 / Ae p=nkF+ n + Ae p=nkF _n_ 



B (k F± ,r,t) = + - P= " KF+ ' + P= ^ F - , (A3) 

mm* m \ n + — n_ n + — n_ J 



where Ae p is defined as, 



2 

The spin-rotation parameter ^sr is obtained by the numerically integration; 

1 r d 2 k 



^SR 



irhnM 2 J (2tt) 



/ vA k2 «■ - "!-) (A4) 



The collision terms are given by the same expressions as given in Ref. |TT| except that the 
quasiparticle energies in the RPA are used. 



8 



REFERENCES 



* E-mail Address: takahasy@eie.yz.yamagata-u.ac.jp 

1 S. Datta and B. Das, Appl. Phys. Lett. 56, 665 (1990). 

2 J. M. Kikkawa and D. D. Awschalom, Nature 397, 139 (1999). 

3 R. Fiederling, M. Keim, G.Reuscher, W. Ossau, G. Schmidt, A. Waag, and L. W. 
Molenkamp, Nature 402, 787 (1999); Y. Ohno, D. K. Young, B. Beschoten, F. Matsukura, 
H. Ohno, and D. D. Awschalom, ibid. 402, 790 (1999). 

4 Ben Yu-Kuang Hu and Karsten Flensberg, Phys. Rev. B 33, 10072 (1996). 

5 K. Hirakawa and H. Sakaki, Phys. Rev. B 33, 8291 (1986); T. Saku. Y. Horikoshi and Y. 
Tokura, Jpn. J. Appl. Phys. 35, 34 (1996). 

6 During the preparation of this manuscript we find that this point is also discussed in Irene 
D'Amico and Giovanni Vignale, Phys. Rev. B 62, 4853 (2000). 

7 A. J. Leggett and M. J. Rice, Phys. Rev. Lett. 20, 586 (1968); A. J. Leggett, J. Phys. C 
3, 448 (1970). 

8 J. W. Jeon and W. J. Mullin, J. Phys. (Paris) 49, 1691 (1988); A. E. Ruckenstein and 
L. P. Levy, Phys. Rev. B 39, 183 (1989); W. J. Mullin and J. W. Jeon, J. Low Temp. 
Phys. 88, 433 (1992). 

9 N. Masuhara, D. Candela, D. O. Edwards, R. F. Hoyt, H. N. Scholz, D. S. Sherrill, and 
R. Combescot, Phys. Rev. Lett. 53, 1168 (1984); J. R. Owers-Bradley, H. Chocholacs, R. 
M. Mueller, Ch. Buchal, M. Kubota, and F. Pobell, ibid. 51, 2120 (1983); L. R. Corruccini, 
D. D. Osheroff, D. M. Lee, and R. C. Richardson, ibid. 27, 650 (1971); J. Low Temp. Phys. 
8, 229 (1972). 

10 L-J. Wei, N. Kalechofsky and D. Candela, Phys. Rev. Lett. 71, 879 (1993); J. H. Ager, 
R. M. Bowley, R. Konig, and J. R. Owers-Bradley, Phys. Rev. B 50, 13062 (1994). 

11 Yutaka Takahashi, Kosuke Shizume, and Naoto Masuhara, Phys. Rev. B 60, 4856 (1999). 

12 In the present study, a typical hydrodynamic scale is few micrometers if we assume optical 
excitations. A microscopic scale is < 100 nm determined by i) the interaction range, ii) 
the interparticle spacings, iii) the mean free path of electrons. 

13 For general references to quantum transport, see L. P. Kadanoff and G. Baym, Quantum 
Statistical Mechanics (Benjamin, New York, 1962); D. C. Langreth and J. W. Wilkins, 
Phys. Rev. B 6, 3189 (1972); H. Haug and A. -P. Jauho, Quantum Kinetics in Transport 
and Optics of Semiconductors (Springer, Berlin, 1996). 

14 G. D. Mahan, Many-Particle Physics 2nd ed. (Plenum, New York, 1990). 

15 H. Haug and S. Koch, Quantum Theory of the Optical and Electronic Properties of Semi- 
conductors (World Scientific, Singapore, 1990). 

16 B. Vinter, Phys. Rev. B 13, 4447 (1976). 

17 S. DasSarma, R. K. Kalia, M. Nakayama, and J. J. Quinn, Phys. Rev. B 19, 6397 (1979). 

18 T. M. Rice, Ann. Phys. 31, 100 (1965). 

19 J. A. McLennan, Introduction to Nonequilibrium Statistical Mechanics (Prentice Hall, New 
Jersey, 1989). 

20 A. E. Meyerovich, Phys. Lett. 107A, 177 (1985). 

21 C. Hodges, H. Smith, and J. W. Wilkins, Phys. Rev. B 4, 302 (1971); A. V. Chaplik, Zh. 
Eksp. Teor. Fiz. 60, 1845 (1971) [Sov. Phys. JETP 33, 997 (1971)]. 



9 



L. Pfeiffer, K. W. West, H. L. Stormer and K. W. Baldwin, Appl. Phys. Lett. 55, 1888 

(1989); V. Umansky, R. de-Picciotto and M. Heiblum, ibid. 71, 683 (1997). 

See ch.29 of N. W. Ashcroft and N. D. Mermin, Solid State Physics (HRW, Philadelphia, 

1976). 

G. Baym and C. Pethick, Landau Fermi-liquid theory (Wiley, New York, 1991). 



10 



FIGURES 

FIG. 1. The inplane energy dispersions of interacting, spin-polarized 2D electrons with a total 
sheet density 2 x 10 11 cm -2 calculated in the RPA. The energy is measured relative to the zero 
of a noninteracting electron dispersion. The Fermi surfaces for spin-up (spin-down) electrons are 
at 0.1127 nm" 1 (0.1115 nm 4 ) for P = 0.1, 0.1373 nm^ 1 (0.0793 nm 4 ) for P = 0.5, and 0.1545 
nm" 1 (0.0354 nm -1 ) for P = 0.9, respectively. 



FIG. 2. The definitions of a longitudinal spin diffusion coefficient D\\ and a transverse spin 
diffusion coefficient D±. See the text for details. 



FIG. 3. The temperature dependence of the longitudinal spin diffusion coefficients -Dy with 0.01 

< P < 0.5. The total electron sheet density is fixed at 2 x 10 11 cm~ 2 . 

FIG. 4. The temperature dependence of the transverse spin diffusion coefficients with 0.01 

< P < 0.5. The total electron sheet density is fixed at 2 x 10 11 cm~ 2 . 

FIG. 5. The magnitude of the effective magnetic field due to the molecular field expressed in 
terms of the spin splitting energy (left) or the Larmor precession frequency (right). 
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